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Abstract 

The charmonium system provides an opportunity to explore a wide variety of 
topics in hadronic physics. Studies of the properties of and transitions among 
cc states yield insight into relativistic and non-perturbative QCD effects. At 
the same time, studies of the decays of charmonium states are a window into 
gluon dynamics and the role of glucball mixing in the production of light quark 
states. A collection of preliminary results utilizing the full CLEO-c ip(2S) data 
sample is presented including two-body branching fractions of Xcj decays, a 
precision measurement of the h c mass, and results on the hindered Ml transi- 
tion: i[)(2S) — > jrjc. 



1 The CLEO-c program 

The CLEO-c physics program is driven by the wide variety of physics accessible 
in the charmonium region. Results from large samples of D and D 8 decays such 



as precision measurements of branching fractions and form factors have direct 
implication on the global heavy-flavor physics program and searches for new 
particles or interactions outside of the Standard Model. The CLEO-c program, 
with the world's largest sample of ip(2S) decays, is also dedicated to exploring 
QCD and hadronic physics through studies of the spectrum and decay of cc 
states below open-charm threshold - this physics topic is the focus of what 
follows. 

The CLEO-c detector, from the beam axis outward, consists of a six-layer 
stereo inner drift chamber, a 47-layer main drift chamber, a ring-imgaging 
Cerenkov detector (RICH), and a CsI(Tl) crystal calorimeter. 1 The entire 
detector is immersed in 1-T solenoidal magnetic field. The two drift chambers 
provide a momentum resolution of ss 0.6% for tracks traversing all layers of the 
chambers. Photon reconstruction with the 7784-crystal calorimeter is achieved 
with an energy resolution of w 5% at 100 MeV and w 2% at 1 GeV. The 
detector covers roughly 93% of the full solid angle. The results that follow are 
derived from a sample of approximately 25 million ip(2S) decays produced at 
rest in the lab by symmetric e+e - collisions in the Cornell Electron Storage 
Ring (CESR). 

2 Two-body XcJ decays 

The three x c j(1P) states are readily produced in electromagnetic (El) tran- 
sitions from the ip(2S) and provide a venue for the study of a wide variety of 
QCD phenomena. Each of the analyses below rely on full reconstruction of the 
entire decay chain: ip(2S) — > 7X c j, Xd — * X. The initial four-momentum of 
the ip(2S) is well known from the beam kinematics; therefore, the experimental 
resolution can be enhanced by doing a four-contraint kinematic fit of the decay 
products to the ip(2S) four-momentum hypothesis. 

2-1 Xc0 : 2 -> 77 

The two-photon decays of the Xc0,2 states are an ideal place to study relativistic 
and radiative corrections to QCD in the charmonium system. To first order 

1 Outside of the magnet flux return is a series of muon chambers; however, 
with a detection threshold of approximately 1.2 GeV/c, the chambers, while 
ideal for B physics, are only marginally useful for charm physics. 




Figure 1: Two-photon invariant mass for \ c candidates. The Xco an d Xc2 peaks 
are visible. The fit to the spectrum (background contribution) is shown by the 
solid (dashed) line. 



these decays are purely QED. Particular interest is in the ratio 1Z = T(xc2 — > 
77) /r(XcO - * 77)i which can be calculated simply at first order as 1Z = 4/15 ^) . 
Deviations from this lowest order calculation result from radiative corrections 
and relativistic effects; therefore, precise experimental determination of this 
ratio is important for validating theoretical calculations that consider these 
effects. 

To measure the two-photon widths, events are selected that contain three 
showers greater than 70 MeV in energy and with | cos#| < 0.75. Events with 
charged particles are vetoed. The signal shapes and efficiencies are derived 
from a signal Monte Carlo (MC) sample that was generated using the nominal 
masses and widths of the \c states ^) . The angular distributions were modeled 
as a pure El transition from the (beam-axis polarized) tp(2S) to the Xc states. 
The decay of the Xc2 state was modeled as a pure helicity-two decay ^). 2 
QED-dominated background shapes are obtained from analyzing data off of 
the ip(2S) resonance at y/s — 3.671 GeV and y/s — 3.772 GeV. 

2 As a check on the systematic uncertainty due to this assumption up to an 
8% helicity-zero component was included, motivated by experimental limits in 
a 2 decay ^) . 



Table 1: Various measured parameters for the decays Xc0,2 — » 77- The er- 
rors are, in order, statististical, systematic, and, where applicable, that due to 
uncertainty in B(ip(2S) — ► 7X00,2) and r tot (x c o,2)- results are preliminary. 

Parameter This Measurement 

B(ip(2S) -> 7Xco) x B( X co -> 77) x 10 5 2.32 ± 0.33 ± 0.15 

B(^(2S) -» 7x c2 ) x S(x C 2 -» 77) x 10 5 2.82 ± 0.29 ± 0.21 

#(x c o -> 77) x 10 4 2.52 ± 0.36 ± 0.16 ± 0.11 

B{Xc2 -» 77) x 10 4 3.20 ± 0.33 ± 0.24 ± 0.18 

T(Xco -» 77) [keV] 2.65 ± 0.38 ± 0.17 ± 0.25 

T(xc2 -> 77) [keV] 0.62 ± 0.07 ± 0.05 ± 0.06 

K = T(xc2 -» 77) /r(xco -> 77) 0.235 ± 0.042 ± 0.005 ± 0.030 



Figure 1 shows the fit to the two-photon invariant mass in data. A total 
of 212 ± 31 and 335 ± 35 events are observed for Xco - » 77 and Xc2 — * 77 re- 
spectively. A summary of the preliminary results appears in Table 1. External 
measurements for B(ip(2S) — > 7Xc./) and r tot (x c j) are needed to obtain two- 
photon branching fractions and partial widths. For these, the PDG average 
is used 2) , and errors due to these external inputs appear as a third, separate 
systematic error in the table. Dominant experimental systematic errors are due 
to the fitting method and determination of the signal efficiency. 

The obtained preliminary result 1Z = 0.24 ± 0.05 is the most precise 
single measurement to date; however, it is not yet precise enough to clearly 
validate any one particular approach for handling both radiative and relativistic 
corrections in the charmonium system. When combined with the PDG value 
of 1Z = 0.18 ± 0.03 one obtains 1Z = 0.20 ± 0.03, a value that only marginally 
disagrees with the zeroth-order prediction of 1Z = 4/15 = 0.27 and motivates 
more careful experimental and theoretical scrutiny. 

2.2 XcJ - V {,) V {,) 

Hadronic decays of Xc states, like J/ip, proceed dominantly through annihila- 
tion into gluons and provide an ideal environment to try to understand gluon 
dynamics and glucball production. Hadronic J/tp decay has received much 
attention in this regard due to the interesting series of results from BES ^) 
that study production of scalar f resonances against flavor-tag ui and <j) states 



Table 2: Measured branching fractions and 90% confidence level upper limits 
for Xc0.2 — > rf'}rf'\ Where applicable, the errors are statistical, systematic, 
and due to B(ip(2S) — > 7Xcj) respectively. All results are preliminary. 



X 


B( X cO ^X)x 10 a 


B( X c2 ^X)x 10 3 




3.18 ±0.13 ±0.18 ±0.16 


0.51 ±0.05 ±0.03 ±0.03 


r/r/' 


< 0.25 


< 0.05 


v'v' 


2.12 ±0.13 ±0.11 ±0.11 


0.06 ±0.03 ±0.004 ±0.004 






(< 0.10) 



in J/ip decay. Close and Zhao interpret the results, which appear to be sug- 
gestive of large OZI rule violating effects, as a signature for scalar glueball 
mixing amongst the /o states in the 1.5 GcV/c 2 region. Following up on this 
work, Zhao proposes a factorization scheme in which one can coherently ana- 
lyze the partial widths of various two-body Xc decays ^) in terms of singly and 
doubly OZI suppressed components, where a large doubly OZI suppressed com- 
ponent could also be indicative of strong glueball mixing. A coherent study 
°f XcJ — * rf^rf'^ provides a testing ground for this production factorization 
model. In addition Thomas notes ^) that these decays provide a mechanism 
to explore the gluonic component of the ?/ . 

Like the two-photon decays, analysis of XcJ — ► rj^rj^ relics on the re- 
construction and kinematic fit of the entire event. This analysis is an update 
of a previous CLEO-c analysis that utilized a subset of the data ^). The 77 
decay candidates are detected in the modes 77, 7r + 7r~7r°, and -fK + ir~, while 
the rj candidates are reconstructed in 77r + 7r~ and r/Tr + 7r~ modes. After the 
kinematic fit, the two body invariant mass distribution is plotted and signals 
are extracted by fitting the peaks to a Breit-Wigner with widths fixed by the 
PDG values ^) convoluted with a MC-determined Gaussian resolution. In cal- 
culating experimental efficiencies, it is assumed that ip(2S) — > 7X c o,2 is a pure 
El transition. Branching fractions are summarized in Table 2. 

The Zhao model casts the branching fractions as a function of r, a pa- 
rameter that is equal to the ratio of the strengths of doubly OZI to singly OZI 
suppressed decays. Figure 2 shows the predicted branching fractions for the 
rjri, r)7]', and 77' 77' final states as a function of r for both the Xco and Xci- One 
notes that for both the Xco an d Xc2 all measured values or limits are consis- 




Figure 2: Predictions for the branching fraction of Xc0.2 to r/rj (solid), r\r\' 
(dashed), and rj'rj' (dotted) as a function the model parameter r, which is the 
ratio of doubly to singly OZI suppressed decays. Experimental measurements 
presented in this work are indicated by the arrows. 



tent with the same value of r lending support for the validity of the model. 
One also notes that the value of r is close to zero which indicates a relatively 
small component of doubly OZI suppressed production, consistent with what 
is commonly accepted as small glueball mixing amongst the isoscalar pseu- 
doscalar mesons. Applications of this technique to the scalar meson sector are 
underway at CLEO-c and BES. 

2.3 Other two-body decays 

In addition to probing the role of glueball mixing in production, two-body 
decays of \ c states provide an opportunity to explore the role of the color octet 
mechanism (COM) in P-wave charmonia. The COM was proposed 

10) 

to 

explain the apparent deficit in theoretically predicated decay rates for nuclcon- 
antincuclcon pairs based on the color singlet model, and therefore motivates a 
new series of precision measurements of two-body P-wave charmonia decay. 

The analysis of these decays proceeds in the same fashion as those men- 
tioned above, namely exploiting the power of the full event kinematic fit to 



Table 3: Measured branching fractions for various two- meson Xcj decays. The 
errors are statistical, systematic, and due to B(ip(2S) — > -JXcj) respectively. 
All results are preliminary. 

~X I B(xco ^X)x 10 a B( X c2 g X) x~IcF~ " 

7T+TT- 6.37 ±0.11 ±0.20 ± 0.32 1.59 ± 0.04 ± 0.06 ± 0.10 

7T 7r° 2.94 ±0.07 ±0.16 ± 0.15 0.68 ± 0.03 ± 0.05 ± 0.04 

K+K- 6.47 ±0.11 ±0.29 ±0.32 1.13 ± 0.03 ± 0.05 ± 0.07 

K° S K° S | 3.49 ±0.01 ±0.15 ±0.17 0.53 ± 0.03 ± 0.02 ± 0.03 



improve resolution and reduce background. The final states reconstructed are 
listed in Tables 3 and 4. The hyperon decays are reconstructed in the following 
modes: A — > pu~ , £+ — » pn , S° — > A7, S~ — > A-7r~, and S° — » A7r°. As 
with the two-photon decays, the kinematically constrained two-body invariant 
mass distributions are fit to extract the yield for each of the Xc signal peaks. 
Experimental efficiency is determined using a MC simulation. For decays of the 
Xd and Xc2 to two hyperons, the helicity of the final state is unknown and the 
range of efficiencies for the allowed helicity configurations is used to quantify 
the systematic error due to this uncertainty. 

The results are summarized in Tables 3 and 4 and represent the most 
precise measurements to date of these two-body branching fractions. While a 
detailed comparison of the results the COM predictions ^) j s no t possible here, 
in general, measured branching fractions tend to be higher than those predicted 
by the COM-motivated predictions suggesting further theoretical understand- 
ing of these decays is needed. 

3 The hindered Ml transition: tp(2S) — ► jrj c 

A clear experimental picture of both the hindered (ip(2S) — > jrj c ) and allowed 
(J/-0 — > -fi] c ) Ml transitions in charmonium is important for understanding a 
variety of theoretical and experimental issues. For example, there is interest 
in using radiative transitions in charmonium to explore photon couplings to 
quarks in lattice QCD Calculating these two rates has been a challenge 
for quark models 

12) _ 

Both rates are key in normalizing exclusive branching 
fractions of the r/ c ; the focus here is on B(ip(2S) — ► jrjc). As will be discussed 
in detail below, the r\ c lincshapc in this decay appears to be non-trivial, and 



Table 4: Measured branching fractions for various two-baryon Xcj decays. The 
errors are statistical, systematic, and due to B(ip(2S) — > -JXcj) respectively. 
All results are preliminary. 



X 



B( X co -^X)x 10 5 B(xd -^X)x 10 5 B( X c2 - X) x 10 5 



PP_ 
AA_ 

s+xr 



25.7 ± 1.5 ± 1.5 ± 1.3 

33.8 ± 3.6 ± 2.3 ± 1.7 
44.1 ±5.6 ±2.5 ±2.2 
32.5 ± 5.7 ± 4.9 ± 1.7 
33.4 ± 7.0 ± 3.2 ± 1.7 
51.4 ±6.0 ±3.8 ±2.6 



9.0 ±0.8 ±0.4 ±0.5 7.7 ±0.8 ±0.4 ±0.5 
11.6 ± 1.8 ±0.7 ±0.7 17.0 ± 2.2 ± 1.1 ± 1.1 

2.1 ±1.4 ±0.2 ±0.1 4.1 ±1.9 ±0.3 ±0.3 
3.3 ±1.8 ±0.2 ±0.2 3.3 ±1.9 ±0.4 ±0.2 

2.5 ±2.1 ±0.2 ±0.2 4.0 ±2.4 ±0.4 ±0.3 

8.6 ±2.2 ±0.6 ±0.5 14.5 ± 1.9 ± 1.0 ± 0.9 



this complicates the measurement of the rate. 

To examine the lineshape in detail, thirteen signal-rich i] c decay modes 
(and charge conjugates) are reconstructed: 2(ir + ir~), 7r + 7r~7r°7r°, 3(tt + it~), 
2(7T+7r-7r ), 2(K + K~), K+K s tt-, K+R-jt , K+R-n+ir-, K + K s k + ti~ k~ , 
K + K^n + 'K^'K a , K + K~2(ir + n~), i]iT + 7r~ , and ■q2(n + 'K^). Like the Xcj decays, 
full event reconstruction and kinematic fitting is employed for these candidates. 
Figure 3 (left) shows the photon spectrum after it has been sharpened by 
the kinematic fit. The background is fit to a (MC-motivated) linear function 
using data in the region > 900 MeV and 560 < E 7 < 600 MeV. Peaking 
backgrounds below 560 MeV are due to h c — ► -fr/ c , photon cascades from ip(2S) 
to Xc to J/tp states, and ip' — > 7r°J/-0, where, for the latter two backgrounds, 
the two photons merge in the calorimeter. The signal shows a distinct tail on 
the high energy side of the photon spectrum. 

Modification of the line shape for this transition is expected since the 
natural width is relatively large and the available phase space grows like E^. 
In addition the hindered Ml transition has an additional Ez term in the matrix 
clement 

13) 

that may enhance the line shape distortion. These additional line 
shape modifications, while theoretically motivated, are not constrained well 
enough to allow a satisfactory fit to the data and lead one to question whether 
this transition is suitable for extracting the mass and width of the r\ c . 

Figure 3 (right) shows the exclusive photon spectrum before kinematic 
fitting (red line) superimposed on the background-subtracted inclusive photon 
spectrum (points). The agreement is excellent indicating that the line shape 
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Figure 3: Left: The photon spectrum the process ?p(2S) —> jr] c for exclu- 
sively reconstructed T] c states and after a full event kinematic fit. The intrinsic 
lineshape dominates the experimental resolution, which varies from 4-7 MeV 
depending on rj c decay mode. Right: The background subtracted inclusive 
photon spectrum (raw shown in inset). The exclusive lineshape before the 
kinematic fit is superimposed as a solid (red) line. 



modification is visible also in the raw inclusive photon spectrum. A variety 
of techniques are used to extract the yield in the inclusive photon spectrum 
including using an empirical parametrization of the peak and simply counting 
events above background. The 10% uncertainty in the number of signal events 
due to uncertainties in background and signal lineshape is the dominant sys- 
tematic uncertainty in the measurement of B(ip(2S) — > 7?7c)- Our preliminary 
result is B(ip(2S) -► 7%) = (4.02 ± 0.11 ± 0.52) x 10~ 3 . 

4 The mass of the h c 

The singlet ft. c ( 1 Pi) state was the last of the expected charmonium states below 
DD threshold to be identified 

14) 

There is interest in understanding the 
hyperfine splittings of the charmonium states as these give insight into the 
nature of the spin-spin interaction in QCD. In the limit that the confinement 
term in the QCD potential carries no spin dependence, one expects non-zero 
hyperfine splitting for only L — states, and the mass of the h c (L = 1) should 



be equal to the spin-averaged XcJ mass. Therefore precision measurement of 
this splitting AM hf (lP) = (M(l 3 Pj)} - M^Pi) provides experimental input 
on the spin dependence of the qq interaction. The error on AMhf is currently 
dominated by error on the mass of the h c ; hence, a more precise measurement 
is desirable. 

The h c is studied in the isospin- violating process ip(2S) — > n h c , h c — > 
7?7 C . Two methods are utilized: one that is inclusive of all rj c decay modes and 
another that reconstructs the r\ c in a collection of exclusive hadronic modes, 
which mostly overlap with those noted in the previous section. Both require 
a signal 7r° from the primary transition ip(2S) — > Tr°h c be identified from two- 
photon candidates within three standard deviations of the 7r° mass and extract 
the signal from fits to recoil mass spectrum against this 7r°. 

The inclusive analysis relies on the identification of a candidate photon 
for the h c — > -yr/c transition with an energy of 503 ± 35 MeV. Removing this 
criteria overwhelms the signal with background and allows one to determine 
the background shape. Figure 4 (left) shows the fitted n° recoil mass spectrum. 
The signal shape is Breit-Wigner with width fixed to that of the Xd convoluted 
with a Gaussian resolution function of width 2.5 MeV/c 2 . The mass obtained 
is 3525.35 ± 0.24 ± 0.21 McV/c 2 . The angular distribution of the photon in 
the h c — ► 7?? c is consistent with that of an El transition: dN/d cos Q 1 a 1 + 
a cos 2 7 , where a = 1. We obtain a = 1.34 ± 0.53 from the data. 

The exclusive analysis, like other exclusive analyses mentioned above, 
relies on full reconstruction and kinematic fit of the entire decay chain. The 
mass of the rj c candidate was required to be within 30 MeV/c 2 of the nominal r] c 
mass ^) . The tt° recoil spectrum from the set of exclusive candidates is shown 
on the right of Figure 4 and is fit to a linear background plus a Breit-Wigner 
convoluted with a double-Gaussian resolution function obtained from MC. The 
mass obtained from the fit is 3525.35 ± 0.27 ± 0.20 MeV/c 2 . 

Accounting for statistical correlations between the exclusive and inclu- 
sive samples, we obtain the preliminary result: M(h c ) — 3525.35 ± 0.19 ± 
0.15 MeV/c 2 . This yields a hypcrfine splitting AM M (1P) = -0.05 ± 0.19 ± 
0.16 MeV/c 2 , which is remarkably consistent with zero. However, Richard 

15) 

cautions against interpreting this result as a lack of spin-spin interactions in 
the IP multiplet as M(h c ) should really be compared with the spin-averaged 
Xc mass as calculated in the potential model, which is several MeV higher. 




Figure 4: Plots of the recoil mass against the 7r° for inclusive (left) and exclusive 
(right) 7] c selection. The h c signal is clearly visible in both cases. Fits to the 
spectra (background contributions) are shown by the solid (dashed) lines. 



5 Summary 

The charmonium system provides a rich landscape to study QCD. We have 
presented precision measurements for many two-body decays of the \cj states, 
which have implications in understanding relativistic and radiative corrections 
in the charmonium system, the role of the color octet mechanism in P-wave 
decay, and glueball mixing amongst the light scalar mesons. The hindered Ml 
transition, ip(2S) — > "/r] c , exhibits a non-trivial lineshape that is necessary to 
understand theoretically if precision experimental measurements for the par- 
tial width for this decay and the mass and width of the rj c are to be obtained. 
Finally we presented a new precision measurement of the mass of the h c that 
is naively consistent with zero hyperfinc splitting in the IP multiplet of char- 
monium. 

Analysis of the large sample of ip(2S) collected with the CLEO-c detec- 
tor continues, and more exciting hadronic physics results are expected in the 
coming year. I would like to acknowledge the work of my CLEO colleagues 
on these analyses. I would also like to thank the Hadron 07 organizers for 
providing a wonderful venue to present these results. 
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